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f^ ■ Abstract. We find two convergent series expansions for Legendre's first incom- 

Cn I plete elliptic integral F{X, k) in terms of recursively computed elementary func- 

5-H ■ tions. Both expansions are valid at every point of the unit square < A, A: < 1. 

Truncated expansions yield asymptotic approximations for i^(A, k) as A and/or 
k tend to unity, including the case when logarithmic singularity A = A; = 1 

fT^ I is approached from any direction. Explicit error bounds are given at every 

order of approximation. For the reader's convenience we present explicit ex- 
pressions for low-order approximations and numerical examples to illustrate 
their accuracy. Our derivation is based on rearrangements of some known dou- 
ble series expansions, hypergeometric summation algorithms and inequalities 
(~| ■ for hypergeometric functions. 
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(<— ^ ' 1. Introduction. Legendre's incomplete elliptic integral (EI) of the first kind is defined by ^ 

\0 : (12.2.7)]: 



F(A,fc)= f —, '^^ (1) 

^ ' J ^(1 - i2)(l - A;2t2) ^ ' 



It is one of the three canonical forms given by Legendre in terms of which all elliptic integrals can 
be expressed. We will only consider the most important case 0<A;<1,0<A<1. 
KA , The subject of series expansions and asymptotic approximations for the first incomplete elliptic 

H \ integral has a long history. An expansion given by E.L. Kaplan in 1948 ^Oj implicitly contained an 
asymptotic approximation for F(A, k) near the singular point A = /c = 1 (see ((3} below). Soon after 
Kaplan's paper various series expansions for the first incomplete EI were given by B. Radon (1950) 
in UH] and R.P. Kelisky (1959) in [TI]- A bit later B.C.Carlson showed in [3] that F{\,k) can be 
expressed in terms of Appell's first hypergeometric series Fi (see |3]), which automatically lead 
to several series expansions through known transformation formulas for F\. In the same paper he 
noted that one can derive rapidly convergent expansions by first expressing Legendre's incomplete 
Els in a different form. This form had later become known as symmetric standard Els. B.C. Carlson 
proved that instead of three Legendre's Els one can use three symmetric standard elliptic integrals 
as canonical forms. The first symmetric standard elliptic integral is defined by [H El El CI : 

oo 

1 /" dt 

Rf(x,v,z) = - / — . (2) 

^ ^ 2 7 ^{t + x){t + y){t + z) 
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It is symmetric in x, y and z, homogenous in all variables of degree —1/2 and related to F{X,k) 

by 

F(A,A;) = Ai^F(l-A^l-A;2A^l). (3) 

Asymptotic formulas for F(A, k) near the point (1, 1) appeared in 1210 El El, but the first com- 
plete asymptotic series with error bounds at each order of approximation was given by B.C. Carlson 
and J.L. Gustafson in terms of the first symmetric standard elliptic integral Rp in |2j. As is clear 
from ^ and homogeneity, the case A, A; — > 1 for F{X, k) is equivalent to the case z ^ oo with 
bounded x and y for RF{x,y,z). The first two approximations from [g] translated into our nota- 
tion read: 



with relative error bound 



41n[(l-fe2A2)/16] 
and 

A 



(2 - X\l + e)) ln(l - k^X'l ^^^^ 2-X\l + k') ^ ^^^ 



n\k) ^ 



+ e2F{X,k), (6) 



with relative error bound 

9{l-k^Xyin{l-k^X^) 3(1 - PX^)^ 

641n[(l-A;2A2)/16] ^^^ 8 ' ^' 

The authors also provide more precise approximations at the price of having the first complete 
elliptic integral in them: 

2 



F(A,t) = -K(Vl-.-).n-^=j^-^=^-*, (8) 



A-(v^r3I^) In _ _ - ^ (2 - A' - .'A' - Vd - A^)(l - P)) + ... (9) 



where absolute errors have bounds given by 

1 - fc^A^ ^ , ^ (l-fc^A2)ln(4) 9(1 - k^X^y 3(1 - fc^A^)^ ln(2) 
^^<^^< Pa^ ' 64 <^^< 2M^ • ^'°^ 

The problem of finding complete asymptotic expansion for z — > cx) with bounded x and y solved by 
Carlson and Gustafson in |6j for the first symmetric standard elliptic integral Rp has been solved for 
the other types of symmetric standard Els by J.L. Lopez in JH]. This research has been continued in 
|14j . where complete asymptotic expansions are found for all symmetric standard elliptic integrals 
when two variables tend to infinity thus settling the question in principle for symmetric standard 
Els. The methods used in the above papers are either based on Mellin transform technique ^ or 
distributional approach J13l I14j . Recent advances in Mellin transform technique can be found in 

As one clearly sees from the error bounds ©, ((TJ and ()1U|) the corresponding approximations 
are only asymptotic when both variables A and k approach one. The approximations derived in 
this paper are of somewhat different character in that only one of the variables needs to approach 
one while the other is allowed to behave arbitrarily including approaching one as well. Hence, 
Theorems El and HI below provide asymptotic approximations for F{X, k) of any order for A and 
k lying on any curve having the endpoint at the side A = 1 or fc = 1 (including the logarithmic 
singularity A = A; = 1) of the unit square [0, 1] x [0, 1] in the (A, fc)-plane. The coefficients of our first 



expansion H42|) are expressed recursively in terms of elementary function. The second expansion (|61() 
contains the first complete elliptic integral minus an elementary function also computed recursively. 
Each approximation is accompanied by a two-sided error bound. Our derivation is based on simple 
rearrangements of certain modifications of some known double series expansions, hyper geometric 
summation algorithms and inequalities for hypergeometric functions (some known and some new). 
The resulting approximations are very precise which is demonstrated in the last section of the paper 
containing numerical examples and a comparison with Q and ©• 

Expansion (|H1|) may be combined with asymptotically precise inequalities for the first complete 
elliptic integral found in j^. These inequalities can be further improved by employing integral 
representations and using the method of refining the Cauchy-Bunyakowsky integral inequality de- 
veloped in [23 - US • Computations show good precision of these results near the singularity. 

2. Expansions of B. Radon and R. Kelisky revisited. In this section we derive two auxiliary 
expansions which will serve as starting points for our main results formulated in sections 3 and 4. 
The first expansion can be viewed as a different guise of a known expansion due to Brigitte Radon 
|18j , while the second follows from an expansion due to Richard Kelisky JJ by some hypergeometric 
transformations. The error bounds found in this section appear to be new. 

To keep the exposition as self-contained as possible we will give an independent derivation of a 
modified Radon's expansion. To this end we need the following lemma. 



Lemma 1 For an integer j > the following identity is true: 

y (l-i^y+i ^ ^ 2j\ '''\l-\)^2j\^^ '> (l-i)„ [l-X'J ' ^''^ 

^-^ 

where the second term is zero for j = 0. Here 

{a)n = a{a + l)(a + 2) • • • (a + n - 1), (a)o = 1, 

is the Pochhammer symbol {or shifted factorial). 

Proof. Euler's integral representation for the Gauss hypergeometric function reads after a 
simple variable change: 

A 


Writing l\12\i for n = j + 1, employing differentiation rule 3, formula 2.1(7)] and representation j3 
formula 2.8(14)] for 2^1(1, 1/2; 3/2; A^) we compute: 



f t^^dt A2J+1 , , , .^ X^^+WdV 





arctanh(-yn) 



u 



=A2 



= r-iyilZ^i /^i + A\ 1 ^^^'^ ^-i)fe+-"(i/2)fe(i/2)^ / A^ y-^- 
( ^v^^/^^^^u(^ + ^^ I ( ^' V ^ Vl'^-^'Vr n-V (i/2).(i/2)n-fc 

^ ^ 2j! Vl-A^^Vl-A^; 2A^^V A2 J ^ ^ l^^ {j - k)k\{n - k)\ 
-( iy^^/2),l /1 + A\ 1 ^^ ..„ (l/2-j)n / A2 y-'^ 



Theorem 1 For A and k satisfying 

1-A:2<(i_a2)/a2, 

and an integer N > the expansion 



(13) 



N 



,^,,,^ = ly^^^^^/'um, 



1 
2A 



j=0 

N j-1 



(j! 



1^2 



(1 - ky 



EE(-i: 



,=ln=0 ii!(l-j)n 

holds true. The bound for the remainder is given by 

A(l/2)7v+i 



(V2wv2-.)„,,_,.,,Y^y-%^,.,.,.) ,H) 



j+n_W_fV2 



1-A2 



\Ri,NiX,k)\ < 



2(iV + l)(iV + l)! 



A2(l-fc2) 
1-A2 



Af+l 



(15) 



Remark 1. It is clear from the error bound (|15|) that expansion (|14j) is asymptotic for (1 
A;)/(l — A) -^ and convergent for A and k satisfying 1)13^ . 



,/2 



U'2+2 



1 + (A:'"t^)/(1 - t^) 



Proof. Put k'^ = I — k"^. Expandin; 
interchanging summation and integration we compute: 

A A 

F{\k)-- 



-1/2 



into the binomial series and 



dt 



i J{l-t^){l-t^ + k'h^) 



dt 



l-t2 



1 + 



k'H-Y'" 



l-t2 



dt /^ , (1/2), k"H^^ 



i:<-)'^*'Vir^rE (-)'^*-/(T 



t^idt 



j=0 -^ i=Ar+i -^ 

Writing the first integral on the right-hand side as Hll|) . we obtain H14|) with i?i,Ar given by 



This series is obviously alternating. The following estimate shows that each term is smaller in 
absolute value than the previous one: 



A 

{j + ir ^ 7(1 



t^^+^dt (1/2), (1/2 + j) 



t2),+2 jl^j + 1) 



(1 - ky 



fj ^2(1 - fc2) 

:i-t2)i+i (i-t2) 



dt 



< 



(1/2), 



j! 



A;2V 



y (1 - t2).+i (1 _ A2) - j! ^' ' W (1 



t^idt 



t2)j+l ' 





The last inequality is due to p3j) . Hence, we are in the position to apply the Leibnitz convergence 
test which implies that the remainder term i?i^7v(A, A:) does not exceed 



t2^+2dt 



(l/2)Ar+i _ J,2-.N+1 

(N + iy.^ ' J (l-t2)iV+2- 



We will prove the following asymptotically exact (as A — > 1) estimate 

A 







valid for all A G (0, 1) and a > 0. Indeed, /i(0) = /2(0) = and 

/{(A) _ 2a 



/^(A) 2a + l-A2 



<1, AG (0,1) 



The estimate ((TB|) immediately leads to (fT^. D 

Remark 2. One can verify that expansion H14|) is a different form of the expansion 






(17) 



due to B. Radon (see HE])- Here Q^- denotes the Legendre function of the second kind. Indeed, 
using representation ()12p and the formula 

Mj + U- + 1/2;; + 3/2; A^) = ^Jzim±lL_Q, (1 

instead of (|lll) in the proof of Theorem ^ we can get (|17l) . 

The function 2-^1 (—f^, 1/2; l;x) creeps up into our considerations on several occasions. It can 
be expressed in terms of Legendre polynomials via 



.F.i-r, 1/2; 1; .) = (1 - .)"/^P„ ^^ - .)^'^ .^(^ - .)'^'^- \^ ^ ,, _ ^.,.,,,^^ (_2_^ 
(see [201 formulas 7.3.1(175)]). Using the first Laplace integral 

P„(z) = - I {z + ^/z"^ -I cos lpY dip 



for P„(z) (see jUl formula (4.8.10)])), we obtain: 

7r 

2Fi(-n,l/2;l;x) =^ f (l - xsin^ ^Y dip. (19) 



We summarize the required knowledge about this function in the following lemma. 

Lemma 2 a) The function Fn{x) = 2-F'i(— ra, 1/2; 1; x) is monotone decreasing for x G [0,1] 
and all non-negative integers n with hounds 

Fn{l) = ^^^ < 2Fi{-n, 1/2; 1; x) < F„(0) = 1. (20) 

n! 

b) For X G [1,2] the function Fn{x) is monotone decreasing when n is odd with bounds 

Fn{2) = < 2Fi{-n, 1/2; 1; x) < F„(l) = ^^^ < 1, (21) 

n! 

and has a single minimum at Xmin £ {^t"^) when n is even with bounds 

Ti 

< 2Fi{-n, 1/2; 1; x) < F„(2) = ^„^^^^^,, < 1. (22) 



c) For X > 2 the function Fn{x) has the sign (—1)" and increasing (decreasing) for even (odd) 
n with the bound 

|2Fi(-n,l/2;l;x)| <(x-ir. (23) 

d) The following identity holds true 

2Fi{-n, 1/2; 1; 1 - x) = ^-^^^Fii-n, 1/2; 1/2 - n; x). (24) 

Proof. There are many ways to prove this lemma. We present a self-contained proof based on 
representation ((T^. 

1. Let < X < 1. Then (UH) shows that F„(x) is decreasing, so that F„(l) < -F„(x) < i^„(0). 
Clearly F„(0) = 1, while F„(l) = (l/2)„/n! is the celebrated Chu-Vandermonde identity and a) 
follows. 

2. Let l<x<2, n = 2A; + l and A; > is an integer. Since 1 — xsin^ ^ is decreasing in x, we 
infer from H19|) that Fn{x) is also decreasing. Hence, Fn{2) < Fn{x) < -Fn(l) and 



F2fc+i(2) = - rcos^''+^^d^ = 



which proves (PT|) . 

3. Let 1 < X <2, n = 2k and A; > is an integer. It is then obvious from (|19() that F"{x) > 0, 
so that F^{x) is increasing. It is also clear that F^{1) < 0. On the other endpoint we have 

F^^{2) = -- rcos^"-V ^~"°"^ # = + - rcos''</'#>0. 
TT Jo 2 vr 7o 

Monotonicity of -^2fc(-^) iii^plies that there is a single minimum at a point Xmin G (1)2) and hence 
< F2k{x) < max(F2fc(l), ^2^(2)). The value of F2k{2) can be computed from (fT^ or found in PUI 
formula 7.3.8(2)]. To prove (|^ we need to show that -F2fc(l) < -^2fc(2) or 

{\)n n\ r(2fc + l/2) r(fc + l/2) 

n! 2"(f!)2 r(2A; + l) r(A: + 1) ' ^ ' 

The last inequality is clearly true for k = 1. Its validity for any integer k follows by induction on 
fc, where the step of induction is secured by the elementary inequality 

(2/c + 3/2)(2A; + l/2) A; + 1/2 
(2fc + 2)(2A: + l) ^ k + l 

Moreover, 



2"(f!)2 r(A; + l 

4. Suppose X > 2, n = 2/c. From (dHl) we see that F2fc(x) > 0, -FaU^) > ^^^ -F2fc(2) > which 
implies that i^2fc(^) ^ ^ ^^^ P2k{x) is increasing. 

5. Suppose X > 2, n = 2A; + 1. As before -F2fc+i(2) = 0, -F2fc+i(^) < ^ and consequently 
F2k+i{x) < and is decreasing. 

6. From the elementary inequality 

|l-ex| < |l-x|, X > 2, 0< e< 1 

and ((T^ we get (PS)) by choosing e = sin^ |. 

7. Finally, identity 1)24(1 is the limiting case of the well-known analytic extension formula for 2F1, 
see formula 15.3.6]. This identity can also be proved by writing F„ as the Legendre polynomial 
as in (UHl) and applying 1201 formulas 7.3.1(175)-(176)]. D 



Theorem 2 For \, k satisfying 
and a positive integer N the following expansion holds true: 



(26) 



F{X,k) = K{k) 



1-A^ 
1-A;2 



1/2 N-l 

E 

m=0 



(1-A 



2\m 



2m + 1 



2Fi(-m, 1/2; 1; (1 - A:")-^) + i?2,iv(A, A;), (27) 



where K{k) is the complete elliptic integral of the first kind. The hound for the remainder term is 
given by 



\R2,N{X,k)\ < 



1 



2N + 1 



{1-X')k 



2\1.2 



-iN 



[1-k^) 



^(l_A2)(l-fc2) 

1 + A;2A2 - 2A;2 



for 1/2 < A;2 < 1, and 



forO< A;2 < 1/2. 



\R2,N{X,k)\ < 



A 



-2 



2iV + l 



:i-A 



2\N 



1-A2 

1-A:2 



(28) 
(29) 



Remark 3. It is clear from the error bound (|2<S|) - (|29|) that expansion (|27|) is asymptotic for 
'1 — A) /(I — k) -^ and convergent for A and k satisfying H26|) . 

Proof. We begin with an expansion around \ = k = given by R.P. Kehsky in JJ| (see also 



i) 



n^ k) = E ^L , ..::, 2F,i-m, 1/2; 1/2 - m; k'^ 



m=0 



(2m + 1) ml 



An application of the identity H24|) transforms it into the expansion 



FiX,k) = Y, 



°o -,2™+! 



A 



m=0 



2m + 1 



,Fi(-m,l/2;l;l-A;2), 



(30) 



valid for fixed < A < 1 and \k\ < 1/A. To make the next step we need the reflection-type relation 

1 



F{X,k) =K{k) 



F(Vl-X^V-ky{l-k^) 



(31) 



which can be easily verified by representing the integral over (0, A) in (^ as the difference of 
integrals over (0, 1) and (A, 1) and introducing the new integration variable u"^ = 1 — t^. The 
branch of \/l — X"^ is chosen so that it is positive for positive values of 1 — A^. The branch choice 
of the second square root is immaterial since F depends on the squared second argument only. 
Expanding the second term on the right-hand side of ()31() into the series ()3Up and splitting the 
resulting series we get (|27j) with the remainder given by 



^2,Af(A,A;) 



1-A;2 



1/2 oo /^ \2\m 

E^^^^^2F,(-m,l/2;l;(l-fc^-- 



m=N 



2m + 1 



To obtain a bound for R2,n we invoke the estimate ()23p. Substituting this estimate into (|32|) yields: 



1-A2 



1-A:2 



|^2,iv(A,A;)| < 
for 1/2 < /c2 < 1 and 

\R2,NiX,k)\ < 

for < A;^ < 1/2. Applying the inequality 

O/m 4- 1 ^ 



1/2 oo 



^^ 2m + l 



m=N 



2^^.2 



(l-A')fc 
(1-A:2) 



l-A^ 
1-A:2 



1/2 oo 



E 

m=N 



(1-A 



2\m 



2m + 1 



m=N 



2m + I 



< 



2S + 2N + 1 - 2N + . 

s=0 s=0 






„Af 



(2A^ + l)(l-a;) 



valid for < X < 1, we arrive at ((2H1) and ^^. U 



3. First asymptotic expansion. Denote 



^^'^^-^'-"m^^-y- 



j=n+l 

The change of summation variable m = j — n — 1 gives: 

^ ^ [(n + l)!]2^ ' ' M 3/2,2 + n,2 + n 



(34) 



We derive a recurrence formula for s„(x) in terms of elementary functions in the following lemma. 

Lemma 3 The functions Sn{x) satisfy the four-term recurrence relation 

4(n + 3)^s„+3 = an{x)sn+2{x) + bn{x)sn+iix) + c„(x)s„(x) + /i„(x), (35) 

where 

an{x) = Sn^ + 36n + 42 - x(2n + 5)^, 

bnix) = 2x{4n'^ + Un + 13) - (2n + 3)^ 

c„(2;) = -4x(n + l)^ 

,^, _ ^(2n + 5)(2n + 3)^ + (n + 3)(8n2 + 24n + 17) ^.^ ,,, ,2. sn+2 
/^ni^J - 8(n + 3)[(n + 2)!]2 ^^''/^^"^ ^ "^^ ' 

and the starting values for the recursion are given by 

so(x) = -21ni±^^±^, (36) 

.,W = (--i)i„ ^ _vTT^ + _ + _ (37) 

^ ^ V 32 4 4y 2 V32 16^ 16 8 64 ^ ' 

Proof. Relation (|35() can be proved by a careful application of Sister Celine's or Zeilberger's 
algorithm ^21- Denote the generic term in ()33() by 



(1/2), (1/2 -i) ,_ J(l/2),f(i-n-l)! 

j!i(l-i)n ^ ^ [j!]2(3/2),_„+i 



The j'-free recurrence relation 

(-Sn^ - 36n - 42)5(n + 2, j + 1) + 4(n + 3)^5(n + 3, j + 1) + x{2n + 5fg{n + 2, j) 

+ 4x(n + l)^5(n, A;) - 22;(4n2 + 14n + 13)5(n + 1, j) + (2n + 3fg{n + 1, j + 1) = (39) 

can be verified by a direct substitution. The difference from the standard algorithms for hyper- 
geometric summation comes from the fact that we have non-standard bounds for the summation 
index (which should be over all integers for standard algorithms). We can, however, remedy this 
by noting that 



Y^ g{n + 2,i + 1) = sn+2{x) -g{n + 2,n + 3), 



j=n+3 

00 



XI 9{n + 3,j + 1) = Sn+3{x), 

j=n+3 



X] 9{n + 2,j) =Sn+2{x), 



i=n+3 



X] 9{n,j) = Sn{x) - g{n,n + 1) - g{n,n + 2), 



j=n+3 



X] 5(^^ + 1,^) = Sn+iC^;) -5'('^ + i,?^ + 2), 



j=n+3 



X 9{n + l,j + 1) = s„+i -5(?^ + l,« + 2) - S'("' + l,n + 3), 

j=n+3 

and summing up the j'-free recurrence ()39p over the range j = n + 3,n + 4, . . .. Together with 
definition of g{n,j) this yields (|35j) . 
To evaluate the initial term 



i=i 



JJ! 



in a closed form note the identities 

1 



1 

X 



VT^^ 



1 



1 

X 



E^-- 



i=o 



E^^--. 



i=i 



J! 



^/r^^ 



(V^.. 



j=l -J-l 



x-' = 2 In 



1 + ^/r^^' 



and substitute x with —x. The expression (|87|) for the next term 



Sim 



(1/2), (1/2 -jI 



. „ (1 



i=2 



(1 - j)n^ 



is derived from the identity 



giV^MV^,. 



E 



.^2 (i-j)ii! ;r^ jji 



(V2), ■ l^ (1/2) 
' 2 £- (j - l)ij! 



^x^ 



and evaluations 



i=2 



3J 



ii/^x^-^ = 21n^ 
Jj! 2 



JJ! 



E 

i=2 
both deduced from (^nj. Hence, 

vJl/^x-^=/flln- 
^(J-I)JJ! 7 W I 



i=i 



^,>-2.A,„^ 



1 + yr^^ 2 



x2 1 + Vl -X 2x 



+ ^/r^ 2t 



l^dt=fl-^Vn 2 



(40) 



(41) 



+ 



xj 1 + VI -2; 1 + VI -2; 



and H41|) follows. Similar but more cumbersome computations lead to formula ()38() for S2(x). D 
The main result of this section is now formulated as follows. 



Theorem 3 For all (A, k) G [0, 1] x [0, 1] and an integer N > 1, the first elliptic integral admits 
the representation 



N ,, ,„x ,, ,„x , N-\ 



Z.M ,^ li 1 + A V- (1/2), (1/2),,, ,2^,- Iv-Vl-A^r /(l-fe2)A2\ ^ .^ ^. 
Fix, k) = -ln^—^Y. j, / (1 - ^')^ + ^ E ( ^A^ ) -« ( S-Ta^ ) + ^^(^' ^)' 

(42) 
where s„(-) is found from H35() - ()38() . T/ie remainder term is negative and satisfies 

^^^^^^^,p /jv+i(A, fc) < -i?7v(A, A:) < 2[(Af + l)!]2 •^^^^' ^^' ^^^^ 

where the positive function 



fNi\k) 



1-q(1-A:2) 



y^- 


1- 


-A2 


+ 1 


qA^( 


1-fc^) 


y^- 


1- 


-A2 


-1 


aA^( 


l-fc2) 



fl-fc2)ln.^ + ^ 



"Av/l + ^A^(TTfc5) v/l + aA2(l-fc2) ~ 1 ) I _ (iV+i/2)2 

(44) 
is bounded on every subset E of the unit square, where 

1 — k 

sup < oo (45) 

k,xeE 1 — A 

and is monotone decreasing in N . 

Remark 4. Error bound (|43|) shows that expansion (|42)) is asymptotic as A: ^ 1 along any 
curve E lying entirely inside the unit square with (|45() satisfied, including those with endpoint 
(1, 1). The expansion is convergent for any fixed 0<A<1,0<A;<1. 

Remark 5. If condition 1)45^ is violated but (1 — /c)'"/(l — A) remains bounded, then m-th and 
higher approximations are still asymptotic. In this case, however, it is much more effective to use 
approximation (|6H) from Theorem |11 

Proof. For the incomplete elliptic integral of the first kind we have 



.„,„.^,„i±^ 5: (!/»,_,., 



^^g|'->-^W^<'-^Kr^r '^^' 



j=l n=0 

according to ()14|) . Rearranging the double sum according to the rule 



oo j — 1 oo oo 

E E ""'J = E E ""'i (^^) 

jr-=Ar+lra=Af n=N j=n+l 



(taken for A^ = 0), we get 



Summing the first series for j from to A^ and the second for n from to iV — 1 and leaving the rest 
as a remainder we obtain (|42jl by Lemma El and definition (|33|) of the functions Sn- The remainder 



10 



term is thus given by 



j=N+l 



2 \ J~" 






n=Af j=n+l 

To estimate Rn we change the order of summations in the second term according to the rule ()47|) 
appUed from right to left. This yields 

j=N+l ^•''' 



1 °° 



j=N+l 



, (1/2), 



1-A2 



i i-i 



(1/2 -i)„ /I -A' 



£^. (l-i)n 



n=N 



2\ ^ 



-A2 



Introducing the new summation variable k = j — n and applying standard hyper geometric summa- 
tion algorithms as realized by Maple "sum" procedure, we get 

\il/2-j)^fl-X'Y (1/2), /i-A2V^^(A:-l)!/ -A^ \^ 



.. (1 



n=N 



(l-i)r 



-A2 



(i - 1)! V -A2 7 ^ (1/2) 



. . (1/2)^ U 



(1/2), 



l-A^ 



2\ i-1 



(j-1)! ("V -A2 

(j-iV)! /1-A2\^-^ 



2i^i l,l;3/2 



k V1-A2 
-A2 



1-A2 



(l/2),_^+i V -A2 



Fi{l,l+j-N;3/2+j-N; 



1-A2 



for the inner sum. Now use 



Fi 1,1; 3/2; 



-A2 



1-A2^ 1 + A 
m 



1-A2y 2A ""1-A 

for the first term in braces and Euler's integral representation 



2F1 l,l+i-iV;3/2 + j-iV; 



A2 \ _ (l/2),_^+i f t^-^{l-t)-y^dt 



1-A2y {j-N)\ J (l + tA2/(l-A2)) 





(l/2),„jv+i A -A 



{j-N)l V -A2 



2\j-N+l 



{-uy-^du 



i (l + ^)A/l-ni# 



for the second. Substituting these expressions into the above formula for Ri\j{\, k) and interchanging 
summation and integration we arrive at 



RNix,k)=t}}i 7 ^\:.^7^1 du ± Ki^[-(i-.2)^F. 



2A J u^{l + u) 




j=N+l 



if-) 



(49) 



An easy computation shows that 



E '^fi-y 



^l/2)Af+lJ f ^N+1 



j=N+l 



(j! 



1^2 



[(iV + l)!]2 



(-x)^'+SF2(1, N + 3/2, N + 3/2; N + 2,N + 2;-x). 
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Hence, 



l-A 



T 



^ ,, ,, (l-A;2)^+i[(l/2);v+i]' f 3F2il, N + 3/2, N + 3/2; N + 2, N + 2; -{I -k^)u)udu 



^^1<" + ')'1^ i (i + „)(i-„-'V'^ 



Next, we apply the inequality 



<3F2{N + 3/2,N + 3/2,l;N + 2,N + 2;-x)< t^ttt^^, (50) 



■^ "•" ~(NTTp~ lNTTp~ 

valid for x > 0. A proof of this inequality will be given elsewhere. Thus we have i?Ar(A, A;) < and 

where 



g{a,X,k) 



a2 



Q [1 + a{l - k^)u] (1 + n) (l - 7xi#) '^' 



1 ' ' lnt^^^^-a-k^)ln'^' 



} (51) 



and ai = [{N + l/2)/(iV + l)]^, 03 = [(iV + 3/2)/{N + 2)]2. Defining 

/^(A, k) = gi[{N + l/2)/(iV + 1)]2, A, A:) 

we obtain the error bound (|15)) . The statement (|^^ about the boundedness of /Ar(A, A:) follows 
from an examination of the right-hand side of ()51() . The monotonicity of /Ar(A, A;) in A'" is implied 
by monotonicity of g{a, A, k) in a which is clear from the integral representation (|5H1 . D 
Remark 6. If simplicity is preferred to precision one can apply the elementary inequality 

< 3^2(1, N + 3/2, N + 3/2; A + 2, iV + 2; -x) < 1, a; > 0, 

instead of H5U|) . Using this inequality and explicit representation 

^^fl-u^) 'udu 2A2 , 1 + A 

Am ■ 



(1 + n) 1-A2 l-A' 





one gets the following error bound: 



< -RMX.k) < "^/^'yj'' -f "" I'-W^' 



[(A + l)!]2 11-A2 2 l-A 

This error bound is relatively precise for k > X but loses precision substantially for X > k, while 
the bound (|^5|) is very precise for all values of parameters. 

The first order approximation obtained from (|42|) (see (|67|) below), 

^,,,, 1,1 + A 1 2 l-fe^l^;^ 

Fi A,A;) = -In r + v In + In -, 

^ ' ' 2 l-A A 1 + ^(1 _ A2A:2)/(i - A2) 8 l-A' 
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has an amazing property to be correct asymptotic approximation for F{X, k) not only as /c — > 1 
but also as A ^ including the case when both A, A: — > along any curve. Indeed one can easily 
check that 



as A — > 0, while 



so that 



m,k) = x+{f^-j-^e + ^k']x' + oix^^ 



F{X,k) = X + l{l + k^)X'^ + 0{X^), A->0, 





F{X,k)-Fi{X,k) = 0{X^), A^O. 



Thus Fi(X,k) is a true approximation for two sides of the unit square (including endpoints) - the 
side A = 0, A; G [0, 1] and the side A; = 1, A G [0, 1]. The same is true for higher order approximations 
but the approximation order for A ^ does not increase with N. 



4. Second asymptotic expansion. Denote 



A„(.)=^i'^^^]J=^iiii^.K 



j=0 



j J{2{n + j) + l)j\ 



(52) 



The following evident formula is more notational than meaningful: 

. / N 1 ^ A 1 3 , 

^■n[x) = s-To —,n-\ — ,n+l;l,nH — ; —x . 



(53) 



We give three representations for An{x) in the following lemma. The first is more convenient for 
computing explicit expressions, the second is designed for easier estimation and the third provides 
a hint for an alternative derivation of Theorem (see details in Remark IH)) . 

Lemma 4 The following identities hold for the functions An{x): 



A.n\x) — D 
n! 



n-l 



(-1)"^ in(yrT^ +^x) + 4^ g(-ir ^#^: 

" ^"~" ■^ (1 — n)j 



,n-j 



n\\ X 







p-1/2 / ^ 

2F1 -n, 1/2; 1; : ) dt 



2nx 



1 



j=0 



P-V2 f 2 + t , , 

Pr. _ r- dt, 



l + tj' 2x"+l/2y (l+t)(n+l)/2 ''\2^/TTt 




(54) 



where the second term, in brackets equals zero for n = and D^ means differentiation in x. 
Proof. From ()52() An{x) can be written in the form 



^ 00 



n] 



1 tmmi,^. 



j=0 



2n + 2i + 1 j! 



(55) 



For a formal power series 



fix) = "^ajX^ 

j=0 



we have {D^ means differentiation in x) 



D:xV{x) = Y.iJ + l)iJ + 2)...iJ + 

j=0 



nja^x- 






(56) 
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and 



1 r °° 

n i=0 



-^^x^-. 



2n + 2j + 1 



(57) 



Putting ttj = {— ly (1/2) j/j\ gives f{x) = 1/^/1 + x by the binomial theorem. Combining (|56|) and 
(|57j) we obtain from H55() : 



E 

i=o 



n + j\ (-1^(1/2), , ^ J_z)",-i/2 /■ ^!:!^ 
i ;(2(n + i) + l)j!'' 2n! ^ 7 ^T+i 



dt. 



(58) 



The integral on the right-hand side can be reduced to pij) by the variable change y^ = t/(l + t) 
t = yy{l-y^),dt = 2ydy/{l-yy: 







dt= I e-\i--^dt = 2 



l + t 



y dy 





n-l 



(1 - y2)n+l 



(59) 



= 2(-i)"ii/^ in( VTT^ + v^) + ^ x:(-i)^'^H^-""^- 

n! nVx ^ (l-^)i 

Hence, from H58() and the above evaluation we arrive at the first formula ()54() . An alternative 
method of evaluating the right hand-side of H58|l is the following. Make the variable change t = ux: 

X 1 

+n-l/2 /■ n-1/2 J 



^/TTI 



Vl + "UX 





differentiate under the integral sign and apply the Leibnitz formula 






D:x-f{x)=nlYr\^D'J{x) 

and the elementary formula 

Z)^(l + ux)-^/^ = (-l)'=(l/2)feu^(l + ux)-''-^/^ 
to get 



2nr^'' J ./TTt 2n\ 








"•yr+ 



nx 



dn 



h. I'"'-'" 



-SOS-'tt^ 



-|- ux 



2n\J 





. fe=o ^ ^ ■ 



ux) 



du 



-fc-l/2 



du 



2 J ^/l + ux 



Finally, substituting back t = ux we obtain: 



I j- n 1/2 , ux , , 

2F1 -n, 1/2; 1; — ) du. 

1 + ux 



An{x) 



1 r t"-V2 / i 

=2^1 -n, 1/2; 1; ^ : ) dt. 



2x"+v2 J ^/r+T 



i + t 







The last equality in (|54|1 is a direct consequence of (|18|1 . D 



(60) 
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Theorem 4 For (A, A;) G [0, 1] x (0, 1] and integer N > 1, the first elliptic integral admits the 
representation 



N-l , 



F{\ k) = K{k) - . _ ' ^ (1 - \yAA-^-^\ + R^{\ k), (61) 

where the functions An{x) are found from 1)54^ . The remainder term is negative and satisfies 

(1 _ A2)^+1/2 



2X^NV2^^W^^ 



> -RN{X,k) > 



{i-xT-'/HmN U ., ,, a-fe^) iJi|oi-A^ v(i-A2)(2-A2-F) 



(62) 



Remark 7. The error bound (|62() shows that the expansion ()61() is asymptotic for A ^ 1 along 
any curve lying entirely inside the unit square, including those with endpoint (1, 1). The expansion 
is convergent for any fixed 0<A<l,0</i;<l. 

Proof. By Theorem |21 for values of A and k satisfying (|26j) we have expansion l|27|) which, after 
little modification, can be written as: 

Changing the order of summation according to the rule 

oo m oo oo oo oo 

/ > / ^ "-mj = 2_^ Z^ 0,i,i-n = ^ 2^ "'n+jjj 
m=0 j=0 n=0 i=n n=0 j=0 

we obtain the formula 

Fix k) - Kik) -Til- X^r T ^' - "^'^'^''^ f " + '^ (-^y^y^^^ (64) 

which in view of H52() . can be split as follows: 

n _ X2U/2 ^~i /I _ \2\ 

F(A, k) = Kik) - ^ _ J ^ (1 - XrAn ( Y3fc2 ) + ^^(^^ fe)^ 

^ ' n=0 ' 

with j4„ defined by (|5^ and 

(1 - A2)V2 



i?iv(A,/c) 






^' ■" ' n=N 



As written in 1)64^. the inner sum does not converge unless k < X. However, it was shown in 
Lemma that Anix) is an elementary function defined for all (A, k) £ [0, 1) x [0, 1). We will prove 
that the outer sum converges for all such A and k. To this end substitute the second formula (|54() 
for An and apply the estimate from above from (|20jl to get: 



-«„(A..) = lE(l-'=Y / ^=^.F,(-„,l/2;l;^)<i. 



n=N 

(l-A2)/(l-fe2) (1_A2)/(1_A;2) 



1 / dt f. /,2uin_ a-fer / t^'dt 

2 y y^(rTt)f^^' 2 y y^(rTo[i - (1 - fe2)t] 
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(1-A2)/(1-A:2) V 2-fc2-A2 

^ (1 - A;2)^ f t^dt _ (1 - A:2)^ /■ y^^dy 



2A2 7 A/t(r+t) A2 J (l-y2)^+i' 



where we used (|59|) to obtain the last equaUty. Now an apphcation of inequaUty (|16j) with a = N 
gives the upper bound in (|H^ . 

To find a lower bound we again apply (|20|) but this time the estimate from below. This yields: 






(1 - A;2)^(l/2)^ f t" 



(l-A2)/(l-fc2) 



2iV! J y^KTTT) 


It follows from .lH Theorem 1.10] that 



/ ^ 2Fl{N + 1/2, 1; AT + 1; (1 - fc^)t)dt. 

J \/t(l + t) 



Consequently, 



2FiiN + 1/2,1- N + l;y)>^^ for yG(0,l). 



(l-A2)/(l-fc2) 
(l_A2)/(l-fc2) 



(1 - fc2)A^(l/2)^ / „ , 



2N\ J V 1 + t 



Both functions t and \ hAri ^^^ non-decreasing for t > 0, A^ > 1. Hence, we are in the position 
to apply the Chebyshev inequality |1H| formula IX(l.l)] which results in: 

(l-A2)/(l-fc') (l-A2)/(l-fe2) 



(l-Ay-V^(l/2)^ /, (1-fc^) 1-A^ V(l-A^)(2-A^-fc^) n 

n 

Remark 8. Expansion (|64|) can be obtained directly from the definition of F{X, k) without a 
use of Theorem 121 Indeed, substituting the last formula ((M|l into (|BHl we get 

1-a2 



F(A,A:)=ir(A;)--y(l-A;2)'^ / ^, ^,,,, Pni Z— ]dt 



n=0 



T3I2- 



'^<^'4 / 7j|fe!S<-'^>"a^^"(^)- '^^* 
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The generating function for the Legendre polynomials is given by |21j : 

^ 1 



Y^Z^Pni^ 



Hence, 



and so 



Vl-2xz + z'^' 

n=0 



(i_fc2)»t" / 2+t \ ^ x/m 



n=0 






F{X,k)=K{k) ^ ' '^^ 



2 J ^t{l + kH){l-t + kH) 



This formula can be obtained from the elementary relation 

1 

du 



Fiyk) = ,m-I 



y/{l-U^){l-k^U^) 
A 



by the variable change t = (1 — u^)/(l — k'^). Thus the whole process could be started from the 
above representation. The change of the integration variable and an application of 1)66^ then give 
the expansion H65() . The representation of the general term of this expansion as n-th derivative as 
in the first formula (|54j) can be then obtained by taking (2 + t)/{2y/l + 1) as a new integration 
variable in (jbSj) and applying the Rodrigues formula for the Legendre polynomials. 

5. Results of computations. In this section we present several examples of computations with 
the expansions obtained above. We also give a comparison with the approximations (JU and (jH)) 
due to Carlson and Gustafson. 

Consider expansion (|42|) first. From (|l-i6|) and (|H7|) we have: 

so ^ 7T^] =21n 



1-A2 J 1 + ^(1-A2A:2)/(1-A2)' 



si 



A2(l-fc2] 

1-A2 



1 



A2(l-fc2) 

2(1 - A2) 



In 



2 1-A:2A2 1 /1-A:2A2 



l + ^(l_A:2A2)/(i-A2) 2(1 -A2) 2V 1 - A2 



Hence, the first and the second order approximations read: 

^/,,^l,l + Al, 2 1-A;2l+A 

FiA,A:) = -ln r + v In + In -, 67 

^ ' 2 1-A A i + ^(i_A2A;2)/(i-A2) 8 1-A' ^ ' 

F,(A,fc) =Fi(fc,A) + fi#-^i-^Vn 



4A 2A3 ; 1 + V(1-A2A;2)/(1-A2) 



il- k^)VT^m? +A(i_,2).i,i±A_ (,3) 



4A^/r^A2 + 4AVl-A:2A2 128' ' 1-A 

Denote by Ajv the difference between the upper and the lower bounds in (|43() : 

Approximation 1)6 7|) combined with inequality ()43() puts F{X,k) within an interval of length Ai, 
while H68() puts -F(A, A;) within an interval of length A2. Numerical results are presented in Table 1. 
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The exact values of F{X, k) shown in the tables below have been computed by Maple wth the 
required number of precise digits guaranteed. 



A 


k 


F{\,k) 


1st order 
approx. (|67|l 


Absolute 
error 


Length of 
error range 

Ai 


2nd order 
approx. (|H5|l 


Absolute 
error 


Length of 
error range 

A2 


.8 


.8 


1.0178 


1.0334 


-.01554 


.742x10-3 


1.0216 


-.00378 


.926x10--* 


.9 


.9 


1.3532 


1.3652 


-.01198 


.657x10-3 


1.3547 


-.00153 


.427x10-4 


.95 


.95 


1.6861 


1.6936 


-.00750 


.430x10-3 


1.6866 


-.4914x10-3 


.143x10-"^ 


.99 


.99 


2.4708 


2.4726 


-.00185 


.107x10-3 


2.4708 


-.2468x10-^ 


.721x10-6 


.95 


.99 


1.7951 


1.7955 


-.405x10-3 


.639x10-5 


1.7951 


-.554x10--^ 


.463x10-^ 


.99 


.999 


2.6240 


2.6240 


-.253x10-* 


.213x10-*^ 


2.6240 


-.350x10-^ 


.157x10-9 



Table 1. Numerical examples for approximations 1^67]) and i|fi8l) obtained from expansion l^42\) . Fifth 
and eighth columns represent differences F{X,k) — Fi{X, k) and F(X,k) — F2{X,k), respectively. The 
numbers Ai, A2 are dehned in If69\) . 

Now we turn to expansion (|61|) . From (|54j) we have 

1 



Ao{x) = —= ln(Vl + X + V^), 



Ax{x) 



1 / 1 . , I ^N \ — X 

— -^ ln(Vl + a; + ^/x) - -j== 



Hence, the first and the second order approximations obtained from H61|) are: 



Fi{X,k) = K{k)-\n\\l 



1-A2 

1-A:2 



1-A2 
1-A;2 



(70) 



l-F 



F2{X,k) = h{\k) - i-^lnf./l + ^4^ + ^^ ^' 



A2-fc2 



l-/fc2 



VI + T^ 



Denote by A^v the difference between the upper and the lower bounds in (|62() : 
(1_A2)^+i/2 (1-A2)^-V2(i/2)^ 



A 



N 



2A2iVV2 - A:2 _ A2 

X W2 - A2 - A:2 



2NN\ 
(1-A;2) 



(71) 



2^/r^^A2 



Jl + 2l^ + 2^-^'''P-^''-^-''' 



1-A;2 



1-A;2 



(72) 



Approximation (|7(H) combined with inequality (|()2j) puts F{X,k) within an interval of length Ai, 
while (jTlf) puts -F(A, A;) within an interval of length A2. Numerical results are presented in Table 2. 



A 


k 


F{X,k) 


1st order 
approx. lf7n|) 


Absolute 
error 


Length of 

error range 

Ai 


2nd order 
approx. ifTTT) 


Absolute 
error 


Length of 
error range 

A2 


.8 


.8 


1.0178 


1.1139 


-.09611 


.1509 


1.0346 


-.01679 


.02932 


.9 


.9 


1.3532 


1.3992 


-.04600 


.0576 


1.3573 


-.00414 


.006075 


.95 


.95 


1.6861 


1.7086 


-.02251 


.0252 


1.6872 


-.00103 


.001387 


.99 


.99 


2.4708 


2.4752 


-.00443 


.0045 


2.4708 


-.408x10-'' 


.5164x10-"^ 


.99 


.95 


2.1496 


2.1523 


-.00271 


.0028 


2.1497 


-.299x10-^^ 


.3102x10-'* 


.999 


.99 


3.0445 


3.0447 


-.200x10-3 


.200x10-3 


3.0445 


-.229x10-*^ 


.226x10-*^ 
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Table 2. Numerical examples for approximations i)70)) and i)7I|) obtained from expansion i|6'I|) . Fiftii 
and eighth columns represent differences F{X, k) — Fi{X, k) and -F(A, A;) — F2 (A, A;) , respectively. The 
numbers Ai, A2 are defined in l[7^) . 

We will compare these results with the corresponding results from "H" given by inequalities © 
and (O. Denote by A*^ and Ag the interval lengthes for absolute error bounds implied by (jSJ and 
O, respectively, i.e. 



A*i = [rhs of © - Ihs of ©]i^(A, k), A*2 = [rhs of - Ihs of ^]F{X, k). 



(73) 



A 


k 


F{\,k) 


1st order 
approx. © 


Absolute 
error 


Length of 
error range 
A*i (see©) 


2nd order 
approx. © 


Absolute 
error 


Length of 
error range 

A*2 (see©) 


.8 


.8 


L0178 


.85814 


.15968 


.2032 


.96415 


.05366 


.12508 


.9 


.9 


L3532 


1.2278 


.12538 


.1304 


1.3291 


.02411 


.05376 


.95 


.95 


1.6861 


1.5993 


.08687 


.0742 


1.6771 


.00900 


.01867 


.99 


.99 


2.4708 


2.4417 


.02910 


.0169 


2.4702 


.647x10-3 


.00115 


.99 


.95 


2.1496 


2.0973 


.05234 


.0409 


2.1466 


.00301 


.00898 


.999 


.99 


3.0445 


3.0306 


.01392 


.0076 


3.0444 


.156x10-3 


.427x10-3 


.95 


.99 


1.7951 


1.7232 


.07182 


.0537 


1.7896 


.00545 


.00750 


.99 


.999 


2.6240 


2.6016 


.02232 


.0115 


2.6236 


.337x10-3 


.368x10-3 



Table 3. Numerical examples for approximations (0) and ^ due to Carlson and Gustafson. Fifth 
and eighth columns equal 9iF{\,k), 02F{X,k), respectively. The numbers A\ A2 are defined in 
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